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ABSTRACT 

In  this  paper  the  eigenwaves  of  periodic  impedance  diffraction  grating  (DG),  reflection 
and  transmittance  on  finite  DG  are  under  investigation.  The  approximate  boundary 
conditions  (ABC)  were  used  for  simulations.  The  refined  expression  for  impedance 
dielectric  strips  is  obtained.  Singular  part  of  approximate  integral  equations  (IE)  is 
extracted  and  analytically  transformed.  As  the  result  the  IE  of  second  type  with  smooth 
kernel  has  been  carried  out  and  solved  with  collocation  method  (CM)  and  Galerkin 
method  (GM). 

INTRODUCTION 

The  application  of  impedance  boundary  conditions  (IBC)  greatly  simplifies  the  solution 
of  boundary  problem  for  structure  consisted  of  fine  dielectric  layers.  It  became  a reason 
of  huge  amount  of  work  appearance  dedicated  as  to  IBC  obtaining,  as  to  their 
application  for  concrete  electrodynamics  problems  solution.  For  example,  the  problem 
of  wave  reflection  from  dielectric  grating  was  solved  in  paper  [1]  with  IBC  method.  The 
main  disadvantage  of  IBC  for  dielectric  structure  is  narrow  frameworks  of  application. 
In  theory  these  conditions  are  correct  if  a layer  thickness  is  rather  smaller  than  a 
wavelength.  In  practice  these  conditions  give  satisfactory  accuracy  if  a layer  thickness  is 
smaller  than  a wavelength.  IBC  application  area  may  be  enlarged,  for  example,  due  to 
modernization  described  in  [2].  In  present  paper  we  elaborate  one  more  type  of  IBC. 
The  IE  described  here  are  obtained  by  means  of  approximate  solution  of  a rigorous  IE 

[3],  Our  IE  become  more  simple  in  the  case  of  E(0,0,E)  wave  diffraction  on  two- 
dimensional  structure  (Fig.  1) 

E(x,y)=Ee{x, y)+k2  jx{x',y')E(x',y')g(x,x',y,y')ds',x,yeS  (1) 

where  EL(x, 0)  is  the  external  field,  x =s-ec,  k - wave  number,  s,ee  - dielectric 
permittivity  of  local  inhomogeneity  and  surrounding  environment,  S - cross-  section  of 
the  local  inhomogeneity,  g(x,x',y,y')  - Green  function  (GF),  in  our  case  the  GF  is  one 
f planar  DW.  The  GF  for  planar  DW  consisted  of  an  arbitrary  number  of  layers  is 
obtained.  Flalf-analytical  solution  for  this  IE  is  already  described  in  [4],  Let's  consider 
the  solution  (1)  for  a single  inhomogeneity.  Process  to  large  quantity  of  inhomogeneities 
is  evident.  Let’s  find  an  approximate  IE  solution  for  rectangular  inhomogeneity. 
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Considering  E{x',y  exp{-ik^£~y'),  where  E(x')  is  unknown  function,  we  have 

]E{x\y')g{x,x\y,y')dy'*  g(xJy.y)]E{x\y')ly'*  g(x, x'y,y)2 

- a -a  ksjE 

And  as  the  result  we  arrive  at  following  one-dimensional  IE. 

/ 

$E(x')3{x,x')dx'  - E(x)/k2 5 = - El’{x,0)/k2b  , 


where  G(x,x')  = ,g(x,x',0,0). 
sin  (k-yfECt) 


5 =2- 


kyfe 


L\  . 


At  k^Ea  « 1 the  expression  for  8 is  well-known 
8 =2  a (e-ec). 


(2) 


(3) 


It  gives  us  an  ability  to  write  the  IE  as  follow 

i 

\j{x’p{x,x')dx'-j{x)/k2 5 = -E1' (x.O)  , (4) 

-/ 

where  ,/(x)  = k28  £’(.v) . We  would  obtain  the  same  IE  if  use  approximate  IBC  from  [4] 
with  supposition  ,/(x)  = -ho  p 0[H(x,a)~  H(x,-a)\.  Considering  8 =oo  in  (4)  yields  the 
IE  for  inhomogeneity  as  metal  strip.  The  current  J(x)  has  a singularity  on  the  border  of 
metal  strip.  One  of  the  methods  to  avoid  this  singularity  is  to  change  variables  as  - 
x = /cos(p,  /((p)  = .//simp  . Numerical  experiments  have  shown  that  the  changing  of 
variables  gives  good  convergence  for  impedance  strips  too,  because  the  current  on  these 
strips  increases  at  the  border  as  well.  Our  IE  (4)  in  this  case  may  be  written  as 

7t 

}/(<p  'M(p^p  'Mp  ' - /(<p  )A 25  / sinq>  =-£■'  (/ cos(p.O) , (5) 

0 

where  g((p,(p')  = G(/cos(p,/cos(p').  Note,  that  the  condition  (p'-xp  yields 
g-(cp,(p')«  g0((p,cp')  = -log|2(coscp -cos(p')j/27r . Then  we  extract  singular  part  of 
kernel  IE  (5)  and  transform  this  IE  as  follow 

71  Tt 

j[/((P’k((P^P')-/((Pko((p,(P,)y(P,+  /t<p)  J^o((P^'Mp'-/(q>)A28  I simp  =-£*■(/  coscp,0)  (6) 

0 0 

The  first  integral  in  (6)  has  no  singularity,  the  second  one  is  equal  to  zero.  Simplified  IE 
(6)  we  solved  with  collocation  method.  To  calculate  the  integral  the  highest  calculation 
accuracy  formula  (formula  of  quadrangles)  has  been  used.  The  second  way  to  solve  (6) 
is  the  Galerkin  method  (GM)  with  cos  /cp ' functions  as  basical  ones.  Upon  transfer  to 
the  x coordinate  they  correspond  to  the  first  kind  of  Chebyshev  polynomials.  There  is  a 
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sense  to  use  cos  m<p  only  for  metal  strips,  for  impedance  ones  the  function  sin  m<p  may 
be  applied. 

RESULTS  AND  CONCLUSIONS 

The  methods  elaborated  have  a good  convergence.  The  number  of  collocation  points 
and  number  of  basical  functions  per  one  strip  is  from  3 up  to  7.  Therefore,  the  wave 
diffraction  by  high  number  of  strips  (up  to  100)  it  is  possible  to  investigate.  The  CM 
codes  calculate  ten  time  faster  then  GM  ones. 

Therefore,  the  eigenwaves  of  periodic  impedance  diffraction  grating  are  investigated. 
Dispersion  characteristics,  the  windows  of  transparency  and  phase  synchronism 
conditions  for  first  and  second  harmonics  (Fig  1 , 2)  have  been  obtained  as  the  result  of 
the  investigation.  Table  illustrates  the  comparison  of  waves  propagation  coefficients 
(for  DG  placed  over  planar  waveguide  h=l;  d/h=0,5;  s =4,0;  £,=4,0;  s2=2,l;  e3 — 1 ,0; 

a/h=0,2;  l/h=0,25;  N=1  is  a number  of  strips)  we  obtained  by  rigorous  method  (variant 

1)  and  by  approximate  one  described  above 
(variants  2 and  3).  Clearly  seen,  that  the 
impedance  we  introduced  gives  higher 

accuracy  for  thick  strips  (ky[za<  1). 
Analogous  results  are  obtained  for 
reflection  coefficient  | Su\  from  dielectric 


Variant 

X 

3,5 

5,0 

6,5 

1 

1,7328 

1,6827 

1,6064 

2 

1,7247 

1,6996 

1,6198 

3 

1,6873 

1,7284 

1,6349 

inhomogeneity  [3], 


ism 
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